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Abstract 

Obtaining initial conditions and parameterizations leading to a model consis- 
tent with available measurements or safety specifications is important for many 
applications. Examples include model (in-)validation, prediction, fault diagnosis, 
and controller design. We present an approach to determine inner- and outer- 
approximations of the set containing all consistent initial conditions/par ameterizations 
for nonlinear continuous-time systems. These approximations are found by occupa- 
tion measures that encode the system dynamics and measurements, and give rise 
to an infinite-dimensional linear program. We exploit the flexibility and linearity of 
the decision problem to incorporate uncertain-but-bounded and pointwise-in-time 
state and output constraints, a feature which was not addressed in previous works. 
The infinite-dimensional linear program is relaxed by a hierarchy of LMI problems 
that provide certificates in case no consistent initial condition/parameterization ex- 
ists. Furthermore, the applied LMI relaxation guarantees that the approximations 
converge (almost uniformly) to the true consistent set. We illustrate the approach 
with a biochemical reaction network involving unknown initial conditions and pa- 
rameters. 



1 Introduction 

The computation of guaranteed inner- and outer-approximations of consistent parameter 
sets of uncertain dynamical systems is important for many applications including model- 
based analysis and verification, system identification, model (in-)validation and controller 
design [151 [H [TBI [H] . 
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We consider the derivation of such approximations for polynomial systems subject to 
unknown-but-bounded (or error-in-variables) data, e. g. measurements. Different methods 
are available in literature to address this problem. For discrete-time systems, for instance 
interval analysis methods, e.g. [H], or relaxation-based methods, e.g. [201 ED El S] can 
be employed. However, both approaches are not directly applicable to continuous-time 
systems without additional assumptions. For instance in [6] it was assumed that the 
time-derivatives of the states are available as measurements, therefore, resulting in a 
steady-state problem similar to [11]. 

One possibility to address continuous-time systems with the mentioned methods is by 
discrete-time approximations, e. g. obtained by numerical integration. However, due to 
the discretization error the consistent parameter sets of continuous-time and discrete- 
time model do not necessarily overlap and, thus, wrong conclusions with respect to model 
validity are possible [22]. A common approach to limit the discretization error relies on 
higher-order Taylor approximations resulting in verified integration methods [191 [9l [TH 
[T3| 121 120], but again the results typically depend on the discretization error. 

A more direct approach uses McCormick relaxations or differential inequalities, see [21] 
and references therein). However, deriving the needed tight state bounds can be difficult. 

Methods using barrier certificates and sum-of-squares (SOS) polynomials [HI [T71 [1] allow 
the continuous-time dynamics to be considered directly without numerical integration. 
However, only few converse results for these, i. e. existence of barrier certificates, are 
known. Furthermore, to the best of our knowledge, no results with respect to approxima- 
tions of consistent parameter sets exist so far. 

The main contribution of this work is the use of occupation measures [H] to derive guar- 
anteed polynomial inner- and outer-approximations of the consistent parameter sets for 
continuous-time nonlinear systems based on results presented in [71[Tni[I2]- The reformu- 
lation in terms of occupation measures leads to a linear but infinite-dimensional decision 
problem. However, its relaxation using truncated moment matrices and their dual SOS 
polynomials is a finite-dimensional convex problem. One particular feature of the em- 
ployed relaxation is the (almost uniform) asymptotic convergence of the approximations 
to the true consistent parameter set. Another advantage is that the continuous-time 
dynamics are completely encoded in the decision problem and, thus, no numerical inte- 
gration is necessary. Furthermore, we exploit the flexibility and linearity of the decision 
problem to incorporate uncertain-but-bounded and pointwise-in-time state and output 
constraints, a feature which was not addressed in previous works. 

This contribution is structured as follows. In Section [2] we formalize the problem setup, 
in particular, the considered system class, the description of the uncertain data and the 
desired properties of the approximations. To obtain constraints for a convex optimization 
problem, we reformulate in Section |3] the polynomial continuous-time dynamics in terms 
of occupation measures. In Section |4| we adapt this approach to the consistent parameter 
estimation problem, which is reformulated as an infinite-dimensional linear programming 
problem. Its solution is approached numerically with a hierarchy of finite-dimensional 
semi-definite programs. We show how to derive outer- and inner- approximations, as well 
as certificates of inconsistency. The approach is demonstrated in Section |5] while its 
advantages and computational issues are discussed in Section |6} 
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2 Problem Formulation 



In this section we state the problem of set-based parameter estimation for continuous-time 
nonhnear systems of the following form 

x{t) = f{t,x{t)), x{to) = xo, (la) 
yit) = h{xit)). (lb) 

Here t G [0, 1] is the time, the states are denoted by x G M""" and the outputs by ?/ G M""^. 
Note that the terminal time is set to one without loss of generality, after a suitable time- 
scaling of the dynamics. Initial conditions (at t = 0) are denoted by Xq. We assume the 
vector field / : M x M"^ — )■ M"^ and h : R x M"^ — > Wy to be polynomial maps. To 
simplify notation, we denote the state vector at time tk G [0, 1] by Xk = x{tk). 

Note that time-invariant parameters p G M"*" can be included in (jl]) by defining states 
with constant dynamics Vz G 1, . . . , n^: 

Xi = fi[t,x{t)) = 0, where Pi := Xi{0) . (2) 

That is, the state vector x contains rix — Up true dynamic variables and Up constant param- 
eters. The parameter values are then given by the initial conditions of the corresponding 
variables Xj. This formulation unifies the tasks of initial condition and parameter estima- 
tion, and also simplifies the notation and analysis using occupation measures and moment 
relaxations (see next sections). 

We assume that constraints on the states, output measurements and initial conditions 
(including the parameters) are given by polynomial inequalities g{-) > 0. We rewrite the 
output constraints gy{y(t)) > as state constraints using the polynomial output map, 
i.e. gy{h(t,x)) > 0. 

Let rrix time-independent constraints on the states be given in the form x{t) G A:',Vt G 
[0, 1] with 

{x: g.^x) >0, V2 = l,...,m,} C M"^ (3) 

We additionally assume that we have a finite set of rrit distinct measurements at time 
points tk, k = 0, . . . ,mt — I, such that: = to < < ' " ' < tmt-i = 1- These time points 
include both the measurements of the output y{tk) and a priori information on the set of 
parameters and initial conditions. For each time point tk, let the constraints on the states 
be given in the form Xk G Xk, k = 0, . . . ,mt — 1 with 

Xk := |x : gxk,i{x) > 0, Vz = 1, . . .,m^^^ C X (4) 

Using this information, we want to estimate the consistent initial conditions and param- 
eters. 

In the following, we denote the admissible state trajectory (an absolutely continuous 
function of time) of ([T]) starting at fixed xq with x{t\xo), i.e. 

x{t\xo) = Xo+ / f{T,x{T)) dr. (5) 
Jo 
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Using ([5]), we define the set of consistent initial conditions based on tlie set of consistent 
trajectories as follows: 



With these notations, we can state the following problem: 
Problem 1 (Consistent parameter estimation) 

Find inner- approximations I and outer- approximations O of consistent initial conditions 
(parameters) such thatX C Xq C O. 

For many applications such as model validation and fault detection, it is sufficient to 
determine if the set Xq is empty or not, we therefore state the following related problem: 

Problem 2 (Certificate of inconsistency) 

// Xq is empty, find a certificate of emptiness. 

Consistency of a single initial condition xq can easily be checked numerically by solving 
equation ([s]). However, to determine the complete set X^ is in general difficult, e.g. 
due to the involved nonlinearities and nonconvexities. We use occupation measures (see 
Section [s]) and convex relaxations to derive the outer-approximation O based on [7J . The 
basic idea to find the inner-approximation X is to determine guaranteed enclosures of 
initial conditions that violate some constraints, cf. Section |4| and |10]. Note that the 




inner-approximation X is given by the complement in X of the union of these enclosures. 

In the next section, we describe how to deal with the continuous-time dynamics and state 
constraints using occupation measures. This enables us to derive convex optimization 
problems that are used to determine the inner and outer- approximations. 

3 Occupation Measures and Liouville's Equation 

The crucial idea we employ in this work is to reformulate the parameter estimation prob- 
lem in terms of occupation measures. This has two main advantages. First, it allows us 
to consider an entire distribution (or measure) of initial conditions and parameter values 
and not just single points. Second, linear relationships encoding the nonlinear dynamics 
(i. e. trajectories) of the system link these initial measures with corresponding measures 
at the intermediate time-points for which measurement data are available. 

In [7] and the references therein, occupation measures were used to estimate the region 
of attraction. Here we first review this approach, before we extend it for parameter 
estimation in Section HI We show that unknown-but-bounded state constraints at inter- 
mediate time points, not considered in [7], can be handled easily. This is achieved by 




(6) 




- 1 
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partitioning the occupation measure w.r.t. a partitioning of the time interval, cf. Sec- 
tion |4j To determine the unknown occupation measures, a convex problem is derived. 
Albeit infinite-dimensional, the convex problem can be solved efficiently by a hierarchy of 
finite- dimensional relaxations [71 [TUl [12] • 



3.1 Occupation measures 

Let M{A) denote the set of finite Borel measures supported on the set A, which can be 
interpreted as elements of the dual space C{A)', i. e. as bounded linear functional acting 
on the set of continuous functions C{A). Let P{A) denote the set of probability measures 
on A, i. e. those measures fi of M{A) which are nonnegative and normalized to ^{A) = 1. 

Now assume that the initial condition xq is not precisely known, but that it can be 
interpreted as a random variable whose distribution is described by a probability measure 
/io G P{X). We define the occupation measure 

fi{A X B) := j J I^xtsit, x{t\xo))fio{dxo) dt 

for all subsets ^ x i3 in the Borel a- algebra of subsets of T x X, where T C M is a time 
interval and x{t\xo) is as in Here, /^(x) is the indicator function of the set A, which 
is equal to one if x E A, and zero otherwise. 

Note that fi G P{T x A") is a probability measure, and that the terminology occupation 
measure is motivated by the observation that the value f^^{dt, B) = fi(T x B) is equal to 
the total time the trajectory spends in the set B C X. In addition, note that encodes the 
system trajectories, in the sense that if v G C°°{T x X; M) is a smooth test function, and 
A^o = ^xo is the Dirac measure at xq, integration of v w.r.t. fi amounts to time integration 
along the system trajectory starting at Xq: 

(f,/i):= / / v{t,x)fi{dt,dx) = / v(t, x{t\xo)) dt. 

The occupation measure fi can be disintegrated as fi{A x B) = Jj^fik{B\t) dt where con- 
ditional fik is a stochastic kernel, in the sense that for every fixed t G T, fik{dx\t) is a 
probability measure on X describing the distribution of the state x at time tk, and for 
every B E X, 1 1— )■ fik{B\t) is a measurable function on T. 

On the one hand, the introduced measures allow us to consider the whole set of initial 
conditions (/io)- On the other hand, it allows a reformulation of the nonlinear dynamics 



(la) as a linear equation which only depends on the occupation measures (//, /ijj.), as 



shown next. 
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3.2 Liouville's equation 

With these notations, for all sufficiently regular test functions v G C^(T x A';]R) and 
r := [0, 1], it holds that 

- ^ [ [UtMtM)M (7) 

X Jx Jt Jx 

The right-hand-side of the above equation can be rewritten 

d \ 
-v{t,x{t\xQ)) + gi&d v{t,x{t\xo))- fit, x{t\xQ))jiio{dxQ) dt 




IT JX 

d 




T JX 



dt 

v{t, x) + grad v{t, x) ■ f{t, x) ) fi{dt, dx). 



dt 



To simplify notation, we introduce the Liouville operator £ : C^(T x X) C(T x X) 
as Cv := If + grad ■ / and its adjoint C : C{T x X)' ^ C^{T x X)' such that 
{Cv,fi) = {v,C'^i) for all v e C^iT x X), i.e. C'fi := - dw{fif). 

With these notations, equation ([T]) can be written concisely as 

{Cv,fi) = {v,6ifii)-{v,6o^o) (8) 

for all V G {T x X) , where 6o and 6i refers to t = and t = 1, respectively. Equivalently, 
in the sense of distributions, we can write 

£V = f^i/Wi - ^o/Uo- (9) 

Equation ^ is Liouville's equation and is also called the continuity equation in statistical 
physics or fluid mechanics. 

Whereas the function x & C{T x X) satisfies the nonlinear ordinary differential equation 
(la), the occupation measure fj. G P(T x X) satisfies Liouville's equation which is a 



linear partial differential equation (PDE) in the space of probability measures. 

Note that as the initial conditions are not known, the measures are unknown as well. 
In the next sections we derive an optimization problem that allows us to determine the 
unknown measures. 

3.3 Estimating the region of attraction 

In [7j it was proved that the region of attraction Xq , defined as the set of initial conditions 



Xo consistent with the dynamics (la) and the constraints x{t) E X , t E T and a constraint 
A*! at t = 1, is the support of the measure /xq solving the infinite-dimensional linear 
programming (LP) problem 

sup (1,/io) 

S.t. /to + /^O = ^1 j^j^g-j 



+ 5oiio - 6ifii = 0, 
Ao > 0, /io > 0, /ii > 0, > 0, 
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where A is the Lebesgue measure restricted to X, i. e. the standard n^-dimensional volume. 



The supremum in (10) is w.r.t. measures flo G P{X), /xq € P{X), /ii G P{X) and 
/i G P(T X X). Note that the slack measure fiQ results from the inequality /xq < A as 
further explained in [7]. The above LP problem has a dual LP 

inf(fo, A) 

s.t. Vo{x) > 0, Vx G X, (11) 

i;o(x) > 1 + t;(0,x), Vx G A", 

vll,x)>0, \/xeXi, 

-Cv{t,x)>0, y{t,x)erxX, 

where the infimum is w.r.t. continuous functions Vq G C{X) and v E C{T x X). 



The above LPs (10) and (11) are infinite-dimensional, because the equations are required 
to hold for all test functions v. One can solve these LPs by a converging hierarchy of finite- 
dimensional linear matrix inequality (LMI) problems using semidefinite programming. At 



a given relaxation order d, the primal LMI is a moment relaxation of primal LP (10) 



whereas the dual LMI is a polynomial sum-of-squares (SOS) restriction of dual LP (11) 
3.4 Sum-of-squares relaxation of the infinite dimensional LP 



The dual LMI w.r.t. (11) is given by 



inf VqI 

s.t. -Cv{t, x) =p{t, x) + qo{t, x)t{l - t) + Y.T=i Qii^^ x)9x,i{^)^ 

Vo{x)=v{0, x) + 1 + ro(x) + Y^'P^^ ro^i{x)gx,i{x), (12) 
Vo{x) =Po{x) + YlT=i Qo,ii^)9x,ii^), 
v{l,x)=pi{x) + Y'^J-^ qi^i{x)gx,^i{x), 

where / is the vector of Lebesgue moments over X indexed in the same basis in which 
the polynomial vo{x) with coefficients Vq is expressed. The minimum is over polynomials 
v{t,x) and Vq{x), and polynomial sum-of-squares p{t,x), qo(t,x), qi(t,x), Po{x), qo^i{x), 
Pi{x), ro(x), ro,i(x), Vi = l,...,m^ and gi^j(x),Vi = l,...,m2.j, of appropriate degrees. 
The constraints that polynomials are sum-of-squares can be written explicitly as LMI con- 
straints, and the objective is linear in the coefficients of the polynomial ^0(2;). Therefore, 



problem ( 12 ) can be formulated as an SDP. 

From the solution of the dual LMI of order d, we obtain a polynomial Vq{x) of given 
degree 2d which is such that the semialgebraic set O'^ := {xq : Vq{x) > 1} is a valid outer- 
approximation of the region of attraction X^, i.e. Xq C O'^. Moreover, the approximation 
converges in the Lebesgue measure, or equivalently almost uniformly, in the sense that 
limrf^oo A(C'^) = X{X*), see [7j for details. 



4 Consistent Parameter Estimation 

We use now the occupation measure approach to address the set-based parameter esti- 
mation Problems [1] and [2j As shown next, this requires several extensions to be able to 
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consider the unknown-but-bounded state constraints at the different measurement time- 
points tfe. 

First, we spht the solution {x{t), t E [0, 1]} of problem ([T]) into nit — l arcs {x{t), t G Tk] 
with Tk '■= [tk,tk+i], k = 0,1, ... ,mt — 2. Now consider their respective occupation 
measures iik,k+i G M{Tk x X), with intermediate measures ^k G P{Xk). Obviously 
Xlfcio ^ = A* ^-iid considering Liouville's equation ^ on each arc of the trajectory 
yields a system of linear PDEs 

jC'nk,k+i = 5t^^^iXk+i - Kl^k, k = 0,l,...,mt~2. 



4.1 Outer- approximation 



An outer-approximation O 3 is given by the support of the measure /io solving the 
LP 

sup (1,/io) 

S.t. flo -\- = X, 

^'fJ'k,k+l = — ^tkf^k, (13) 

Ao > 0, /io > 0, 

yUfc+i > 0, /ifc,fc+i > 0, k = 0,1, ... ,mt - 2, 

where the supremum is w.r.t. measures fiQ G P{X), jiQ G P{X), ^k+i £ P{.'^k+i)^ 
/^fc,fc+i £ ^([0, 1] ^ X), k = 0,1, ... ,mt — 2. The above LP problem has a dual LP 



inf {vq,\) 
s.t. fo(x) > 0, 

Vo{x) > 1 + Vq^i{0,x), 
Vk-l,k{tk,x) >Vk,k+l{tk,x), 
Vmt-2,mt-l(tmt~l^x) ^ 0, 



Vx G A', 
Vx G A", 

Vx G A'fc, /c = 1, . . . , — 2, 

Vx G Xmt-l, 



(14) 



-£vfc,fc+i(t, x) > 0, 



V (t, x) G r X A", A; = 0, 1, . . . , - 2, 



where the infimum is w.r.t. continuous functions Vq G C{X), Vk.k+i £ C'([0, 1], '^), k 
0,1,..., mt- 2. 



As in Section 3.3 , the above infinite-dimensional LPs ( 13 ) and ( 14 ) are solved by a converg- 
ing hierarchy of finite-dimensional LMI problems. From the solution of the dual LMI cor- 
responding to ( 14 ) we obtain again a polynomial Vq{x) such that O'^ := {xq '■ Vq > 1} ^ Xq 
and \imd^^\{xi) = \{X*). 



4.2 Inner-approximation 

For the inner- approximation X C Xq, we build on an idea suggested in [10]. However, we 
have to take care of the different measurements at time-points tk. 

In the following, we consider the set of initial conditions Cq for which there exists an 
admissible trajectory ([s]) that violates at least one of the constraints that define X and 
Xk. By continuity of solutions (the vector field / is polynomial and Lipschitz on the 
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compact set X), this set is equal to 



Xq : 3x(t\xo) s. t. 

._ I 3t G [0, 1] and 3iy s. t. gx,v{x) < 

° ' ' or 

3^^ and 3?7 s. t. gx^^nix) < 



(15) 



where u = 1, . . . , rrix and 77 = 1, ... , m^;^ index the violated constraint, and k = 0, . . . , rrit — 
1 describes the time-points. Obviously Xq := X \Cq. 

Note that depending on x{t\xo), there are different combinations of constraints (gx^,r] 
and gx,u) that can be violated. We directly deal with the different combinations and 
derive for each one an outer-approximation of the set of initial conditions (and hence 
parameters) that lead to the violation of the constraint. As will be detailed in the sequel, 
the inner-approximation is then obtained from the union of complements of these outer- 
approximations. 

To simplify the presentation, we assume that the constraints defining X are not violated, 
i.e. x(t\xo) G X,\/xo G Xo,\/t G T. This is a mild assumption, since the bounds X 
can often be derived from system insight (e. g. mass conservation in chemical reaction 
networks), or can be chosen sufficiently conservative. In any case, the constraints defining 
X can be treated similarly to Xj^. 

Note that the number of possible combinations, rric = 2'"*'™"='= — 1, can be very large. 
We can reduce the number of combinations significantly due to the observation that if 
one constraint at is violated, then it does not matter if the constraints for A; > k are 
satisfied or not and can therefore be ignored. This can be formalized by: 

Cn,ri ■= (16) 



'^k = '^k, \/ k < K 

x^ = h- 1, k=.. 

( gx.,r,{x) < J 
Xk = X, yk>K 

where i = 1, . . . , rrix^. 

Remark 1 (Strict and non-strict inequalities) In equation |7^ , we consider strict in- 
equalities. To account for strict inequalities small numbers (slack variables) are introduced 
when the LMIs are solved. 

Then, once the m^^j. -mj different outer-approximations 0{Ci^^ri) have been determined, an 
inner-approximation is obtained by 

X:=X\ U 0(C.,,). (17) 

K=l,...,mf — 1 
7?=l,...,m^, 
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4.3 Inconsistency certificate 



Solving Problem[2| i. e. certifying emptiness of the set of consistent parameters X^, was not 
addressed in [7j. Mathematically, this amounts to certifying infeasibility of the infinite- 
dimensional LP 



find /to 5 /^fc 5 /^mt -1 5 f^k,k+l 

s.t. /to + /io = A, 

Ao > 0, (Uq > 0, 

fJ'k+1 > 0, Hk,k+l > 0, 



0,1,. 



which corresponds to (13) without a cost function. We can check that we meet all the 



assumptions to apply the generalized Farkas theorem of [5l Theorem 2] and that non- 
existence of measures fio, fJ'k, f^mt-iy f^k,k+i solving LP problem (18) is equivalent to the 
existence of continuous functions Vo,v solving the dual LP 



find ftQ, flk, fJ-mt-li f^k,k+l,k — I, 

s.t. {vq, A) = -1, 
Vo{x) > 0, 

Vo{x) > 1 + t;o,i(0,x), 

Vk-l,k(tk,x) > Vk,k+l(tk,x), 
Vmt-2,mt-l{tmt-lj ^ 



Vx G Xk, k = 1, 



Vx e X, 



(19) 



[tmt-l,X) > 0, V X 

-Cvk,k+i{t, x)>0, W{t,x) eTx X,k = 0,l,... 



nit 

Vx e X, 



■nit — li 

,mt - 2, 



If Xq is empty, then LP problem (19) is infeasible. If an LMI relaxation of problem 



( 19 ) is infeasible at some order d, certified by a normalized Farkas vector in the dual LMI 



relaxation (cf. Section 3.4), LP is also infeasible. Thus, finding a normalized Farkas vector 



thus gives a sufficient condition that can be used to provide a certificate of inconsistency. 



5 Example 

We consider a biochemical reaction network in which the substrate xi is enzymatically 
converted into a product via intermediary complex X2 [25] . The continuous-time dynamics 
is given by: 

X2{t) = PlXi{t){l - X2{t)) - {p2 + Ps) X2{t). 



The following constraints on the parameters and states were used: 

X : 



(pi - 0.1){10 - Pi) > 0, z = l,2,3, 
Xi(l-Xi) > 0, z = l,2. 

Measurement constraints Xk were generated from a simulated nominal trajectory with 
Pi = P2 = P3 = 5.05, xi(0) = 0.9 and X2(0) = 0.05. An absolute uncertainty of ±0.025 
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Figure 1: Set-based estimation of consistent parameters and initial conditions. Consistent 
samples (black dots) were determined by uniform sampling from X and subsequent nu- 
merical integration. Converging hierarchy of outer-approximations O"^ with ci = 3,5,8,9 
for estimation of a;i(0), X2(0). 



was added to the nominal values at the measurement time points to = 0,ti = 0.3,^2 = 1 
to obtain A'o, X2. Thus, at each time-point four constraints were used, which results 
in nic = 4-3. To avoid numerical troubles, note that the dynamics (20) and X were scaled 
such that the values of the parameters range from to 1. 

As can be seen in Figure HI the sequence of outer-approximations O converges to the 
consistent parameter set for increasing d. Using YALMIP and Sedumi, the solving time 
was about 5 seconds for 0? = 3, and 40 minutes for d = 9. 

The inner-approximation in Figure [2] corresponds to the complement (in A") of the union of 
twelve outer-approximations 0{Ci). However, only eight outer-approximations are shown 
since the other four were empty. Solving time was on average 5 minutes per problem. 

In addition we determined the outer-approximation of consistent initial conditions a;i(0), ^2(0) 
and parameters pi,p2,P3, see Figure [3| Solving time was about 4 hours. 



Note that the inconsistency certificates derived in Section 4.3 could also be used to inval- 
idate entire regions in the space of the parameters and initial conditions. 



6 Discussion and Conclusions 

Occupation measures are a classical tool in kinetic theory, statistical physics, optimal 
transport, Markov decision processes, amongst others. In a broad perspective, the po- 
tential of occupation measures and subsequent convex relaxations are not often used in 
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0.15 




Figure 2: Set-based estimation of consistent parameters and initial conditions. Consis- 
tent samples (black dots) were determined by uniform sampling from X and subsequent 
numerical integration. Inner-approximations (white area) results from the union of outer- 
approximations (grey area and black lines) for d = 7. 



systems control. We used occupation measures to approximate consistent parameter sets 
for continuous-time nonlinear systems without the need of numerical integration. A par- 
ticular feature of the derived approximations is the (almost uniform) convergence to the 
true (possibly nonconvex) consistent parameter set. As demonstrated at the example, 
outer- and inner-approximations can be obtained even though only few measurements 
with relatively large error are used. Tighter approximations are expected if more mea- 
surements are used, or if e. g. outer-approximations are iteratively used to refine the results 
as it proved useful for linear relaxations [26| |2T| [3] . 

While inconsistency certificates can be used to prove inconsistency of entire models or 
parameter regions, the outer- and inner-approximations can be used to get the shape of 
the consistent parameter set. Such a description of the outer-approximation by a single 
polynomial can be very useful in some applications. In other cases different types of 
representations like a collection of half-spaces might be more beneficial. The inconsistency 
certificates could also be used to derive in an iterative and recursive manner either outer 
bounding boxes or a description of the consistent set using a bisection algorithms (cf. 
[26]). 

A drawback of the presented approach is the computational workload resulting from the 
LMI constraints. Theoretically, the resulting problems can be solved in polynomial time 
w.r.t. the input size. However, due to the size limitations of state-of-the-art SDP solvers, 
this approach is at the moment restricted to problems of small dimensions. An alternative 
could be the use of LP relaxations for larger dimensional systems. Note that the specific 
geometry of the LMI constraints makes these relaxations typically much more accurate 
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Figure 3: Set-based estimation of consistent parameters and initial conditions. Consistent 
samples (black dots) were determined by uniform sampling from X and subsequent nu- 
merical integration. Outer-approximations id = 3) of consistent initial conditions a;i(0), 
X2(0) and parameters Pi,P2,P3 projected numerically onto the parameter space. 



than the LP relaxations and a trade-off between accuracy and problem size has to be 
made, see in particular the discussion in [121 Section 5.4.2]. 

As an interesting extension, one could consider statistical information, i. e. further con- 
straints on the moments of the occupation measures. Here we assumed no more infor- 
mation such as statistics or probability distributions to be given. In many applications 
where there is a limited number of replicates, i. e. too few to obtain a meaningful statistic, 
this is actually the case. However, statistical information can be included if the data are 
polynomial or information about moments are available jT2| [23]. 
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